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Abstract. This paper concerns a diffusive logistic equation with a free boundary 
and sign-changing intrinsic growth rate in heterogeneous time-periodic environment, in 
which the variable intrinsic growth rate may be “very negative” in a “suitable large 
region” (see conditions (HI), (H2), dO])). Such a model can be used to describe 
the spreading of a new or invasive species, with the free boundary representing the 
expanding front. In the case of higher space dimensions with radial symmetry and the 
intrinsic growth rate has a positive lower bound, this problem has been studied by Du, 
Guo & Peng m- They established a spreading-vanishing dichotomy, the sharp criteria 
for spreading and vanishing and estimate of the asymptotic spreading speed. In the 
present paper, we show that the above results are retained for our problem. 
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1 Introduction 

In the real world, the invasion of harmful species and/or introduction of new and beneficial species 
are natural phenomena. It is of primary importance to predict and analyze the growing and 
spreading mechanism of biological invasions. A lot of mathematicians have made efforts to develop 
various invasion models and investigated them from a viewpoint of mathematical ecology. Most 
theoretical approaches are based on or start with single-species models. In consideration of the 
heterogeneous environment, the following problem 

ut — dAu = a{t, x)u — b{t, x)u'^, t > 0, x G D, 

< B[u] =0, t > 0, X G dD, 

_ u(0, x) = uo{x), X G D 

is a typical one to describe the spread, persistence and extinction of the new or invasive species 
and has received an astonishing amount of attention, please refer to m-m and [m Ea [23 [251 EH 
for example. In this model, u{t, x) represents the population density; constant d > 0 denotes the 
diffusion (dispersal) rate; a(t, x) and b{t, x) represent the intrinsic growth rate and self-limitation 
coefficient of the species, respectively; D is a bounded domain of M^; the boundary operator 
B[u] = + a and /? are non-negative functions and satisfy a + fi > Q, v is the outward unit 
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normal vector of the boundary d^l. The corresponding systems with heterogeneous environment 
have also been studied extensively, please refer to [6113 ED nasi] and the references cited therein. 

In most spreading processes in the natural world, a spreading front can be observed. When a 
new or invasive species initially occupies a region Qq with density uo{x), as time t increases, it is 
natural to expect that flo will evolve into an expanding region 17(t) with an expanding front 
inside which the initial function uq(x) will evolve into a positive function u{t,x) governed by a 
suitable diffusive equation, with u{t,x) vanishing on the moving boundary 

In the natural world, for most animals and plants, their birth and death rates will change with 
the seasons, so the intrinsic growth rate a{t, x) and then the self-limitation coefficient 6 (t, x) should 
be time-periodic functions. Especially, in the winter of severe cold and cold zones, animals cannot 
capture enough food to feed upon and do not breed, seeds cannot germinate and buds cannot grow 
above ground, so their birth rates are zero. In the meantime, their death rates will be greater. 
Therefore, in some periods and some areas, the intrinsic growth rate a(t, x) may be negative. In 
order to simplify the mathematics, in this paper we only consider the one dimensional case, i.e., 
= 1, and assume that the left boundary is fixed: x = 0. For a more realistic description of the 
growth mechanism and spreading of a new or invasive species, throughout this paper we assume 

(H) Functions a, 6 G (( 7 ^/ 2 - n L°°)([0,oo) x [0,oo)) for some u G (0,1), and are 
T-periodic in time t for some T > 0. The function a{t,x) is positive somewhere in 
[0, T] X [0, oo), the function b(t,x) satisfies bi < b{t,x) < 62 in [0, 00) x [0,oo) for some 
positive constants 61 , 62 - 

Motivated by the above natural phenomena, in this paper we shall study the following free 
boundary problem 


ut — duxx = a{t, x)u — b{t, x)v?, 
0 ) = 0 , u(t,/i(t)) = 0 , 
h'{t) = -nux{t,h{t)), 

_/i( 0 ) =/lo, u{0,x) = uo{x), 


t > 0, 0 < X < h{t), 
t > 0 , 
t > 0 , 

0 < X < ho, 


( 1 . 1 ) 


where, ho denotes the size of initial habitat, g, is the ratio of expanding speed of the free boundary 
and population gradient at expanding front, it can also be considered as the “moving parameter” 
of the free boundary, B[u] = au — fiux, a and [3 are non-negative constants and satisfy a + j3 = 1, 
the initial function uo{x) satisfies 

Mo G <^^([ 0 , ho]) , Mo > 0 in ( 0 , ho), S[mo]( 0 ) = uo{ho) = 0 . 

The free boundary condition h'(t) = —fiUx{t, h{t)) in (II.ip is a one-phase Stefan condition and can 
refer to [3] for the more detailed background. 

When both functions a and b have positive lower and upper bounds, i.e., there exist positive 
constants Ki, K 2 such that 

Ki < a{t,x), b{t,x) < K 2 , 

recently, Du, Guo & Peng in m have studied the problem ( 11 . 11 ) for the higher dimension and 
radially symmetric case. By developing the earlier techniques and introducing new ideas and 
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methods, they have obtained various interesting results. When 6 = 1 and a = a{x) changes sign, 
the problem (II.ip has been studied by Zhou & Xiao [JO] and Wang [36]. Complete conclusions 
about spreading-vanishing dichotomy, sharp criteria for spreading and vanishing and asymptotic 
spreading speed of the free boundary were obtained by the author. 

The main aim of this paper is to study the dynamics of (jl.lj] and show that most results of 
HU [Ml HQ] continue to hold in the more realistic situation that the T-periodic intrinsic growth 
rate a{t,x) may be “very negative” in the sense that both \{a{t,x) > 0}| <C \{a{t,x) < 0}| and 

a{t, x)dx = —oo are allowed for any t > 0. To our best knowledge, the present paper seems to 
be the first attempt to consider the sign-changing intrinsic growth rate with the time-periodic and 
spatial heterogeneous environment in the moving domain problem. 

In the special case that the functions a and 6 are independent of time t and have positive lower 
and upper bounds, such kind of problems has been studied systematically. For example, when a 
and 6 are positive constants, the problem (|l.ll) was investigated earlier by m for a = 0 and by [23] 
for 13 = 0] when a = a{x) and 6 = 6(x), the problem (II.ip was discussed by Du, Guo &: Liang [9l[T2] 
for the higher dimension and radially symmetric case; the non-radial case in higher dimensions was 
treated by Du & Guo m- Instead of u{a — bu) by a general function /(tt), Du, Matsuzawa & Zhou 
m, Kaneko |22] and Du & Lou m investigated the corresponding free boundary problems. 

Peng &: Zhao [33| studied a free boundary problem of the diffusive logistic model with seasonal 
succession. They considered that the species does not migrate and stays in a hibernating status 
in the bad season. The evolution of the species obeys the Malthusian equation ut = —5u in the 
bad season, and obeys the diffusive logistic equation with positive constant coefficients in the good 
season. The diffusive competition system with positive constant coefficients and a free boundary 
has been studied by Guo &: Wu [T8|, Du & Lin m and Wang & Zhao [38|- The diffusive prey- 
predator model with positive constant coefficients and free boundaries has been studied by Wang 
& Zhao [MlEZllM]. 

This paper is organized as follows. In Section 2, we prove the global existence, uniqueness, reg¬ 
ularity and estimate of {u,h). Especially, the uniform estimates of ||it(L OllcboXfi)] for t > 1 and 
\\h'\\Q(l+l,)/ 2 (^^J^_^_l^n+ 3 ]) for n > 0 are obtained directly regardless of the size of hoo ■= limt^oo^(i)- 
These uniform estimates allow us to assert h'{t) —>■ 0 when hoo < oo and play a key role for deter¬ 
mining the vanishing phenomenon. In Section 3, we exhibit some fundamental results, including 
the comparison principle for the moving domain, and some properties of the principal eigenvalue of 
a T-periodic eigenvalue problem. In Section 4, we shall derive the spreading-vanishing dichotomy: 

Either 

(i) Spreading: limi_,.oo= oo and limn^oo uit + nT, x) = U{t,x) uniformly on [0,T] x [0,L] 
for any L > 0, where U{t,x) is the unique T-periodic positive solution to (|4.2p : 

or 

(ii) Vanishing: limt^oo h{t) = hoo < oo and limi_,.oo maxg< 3 ,<;j(j) u{t, x) = 0. 

In Section 5, we first establish the sharp criteria for spreading and vanishing, then estimate the 
asymptotic spreading speed of the free boundary. The last section is a brief discussion. 
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We should remark that for the higher dimensional and radially symmetric case of dni), the 
methods in this paper are still valid and the corresponding results can be retained. 

2 Global existence, uniqueness, regularity and estimate of solu¬ 
tion 

In this section, we give the existence, uniqueness and estimate of the solution {u,h) to (jl.ljl . 

Theorem 2.1 The problem (II.ip has a unique global solution {u,h) and 

/iGC7i+^((0,oo)), (2.1) 

where 

Doo = ■■ t> 0, X e [0,/i(t)]}. 

Moreover, there exists a positive constant M = M (||a, 6, uqIIcxd) such that 

0 < u{t,x) < M, 0 < h'{t) < fiM, Vt>0, 0<x< h{t). 

Further more, if h^o '■= h{t) < oo, then there exists a positive constant C 

such that 

\W\\c‘'/'^([n+l,n+S\) — G V n > 0, (2.3) 

\\u{t,-)\\w([0,h{t)]) <C, V t > 1. (2.4) 

Proof. Noting that the functions a{t, x) and b{t, x) are bounded, applying the methods used 
in Bm with some modifications, we can prove that dnn has a unique global solution {u,h), 
u E (Doo), h E (0, oo) and u satisfies the first estimate of (12.2p . The second estimate 

of ()2.2p can be proved by the similar way to that of [361 Theorem 2.1]. The details are omitted 
here. 

Now, let us prove (12.11) . Let y = x/h{t) and w{t,y) = u{t,x). A simple calculation gives 

Wt - dC{t)wyy - ^{t, y)wy = a{t, h{t)y)w - b{t, h{t)y)uP‘ , f > 0, 0 <y <1, 

' (*> 0) = 0, w{t,l)=0, t>0, (2.5) 

, 'w(0, y) = uo{hoy), 0 < y < 1, 

where C{t) = h~‘^{t), f,{t,y) = yh'(t)/h(t). This is an initial-boundary value problem with fixed 
boundary. Remember a, b ^ (^^’^([OjOo) x [0,oo)), u E and h E oo). For 

any given r > 0 and 0 < e <C 1, applying Theorem 10.1 of [2H Chap.4, p.351] to the problem (|2.5p 
in [s,t] X [e, 1] and [s,t] x [0,1 — e], respectively, we obtain that 

w X [£,l])nC'i+§’2+^([e,r] x [0,1-e]). 

This implies 

ueC^+^’^+''{[£,T] X [£h{t),h{t)]) n X [0,(1-e)h(f)]). 


( 2 . 2 ) 

= C{gL,hoo,M) 



A diffusive logistic equation in time-periodic environment 


5 


Due to the arbitrariness of e, one achieves 

where Dr = {{t,x) : 0<t<T, 0<x< h{t)}. Hence, by the condition h'{t) = h{t)), it is 

immediately to get h' E r]). 

Finally, we prove ()2.3p and (|2.4I1 . For the integer re > 0, let w'^{t,y) = w{t + n,y). Then rc” 
satisfies 

' < - dCw^y - = r{t, y), 0 < t < 3, 0 < y < 1, 

< (^aw^ — ■^Wy'^ (t, 0) = 0, 1) = 0, 0 < t < 3, 

, y) = re(re, h{n)y), 0 < y < 1, 

where C"' = ({i + n), + n, y), = h{t + re), and 

P{t, y) = a{t + re, h'^{t)y)w'^{t, y) - b{t + re, h^{t)y){w'^{t, y))^. 

Noticing (12.2h . we have that re”', and /” are bounded uniformly on re, and 

cu”(r) = max |C”(s) — C"'(t)| < —)• 0 as r —)• 0 

0<s,t<3, |s-t|<r hg 

uniformly on re. Moreover, C”(^) ^ for all re > 0 and 0 < t < 3 as h{t) < hoc < oo. 

Remember the boundary conditions of w^{t,y) at y = 0,1 and the fact 0 < < y,M. 

Choose y S> 1, we can apply the interior estimate to derive that there exists a positive constant 
C independent of re such that, for all re > 0, 

(i) when /3 = 0, we have gj^jQ < C (see [28l Theorem 7.15]); 

(ii) when/3 > 0, we have < C (see [28l Theorem 7.15]), 3 ]x[o t]) 

< C (see [28l Theorem 7.20]). 

In a word, < C for all re > 0. In view of the embedding theorem, it follows that 

llre^”'ll i+i' 1 J. < C for all re > 0. This implies llrcll i+v , , < C for all re > 0, where 

" "c^’^+‘'([l,3]x[0,l]) “ “ n II 11 ^^,_ _ , 

En = [re + 1, re + 3] x [0,1]. This fact combined with 

h'{t) = Ux{t,h{t)) = h~^{t)wy{t,l), 0 < h'{t) < yM, 

allows us to derive (IZ3I). Since these rectangles En overlap and C is independent of re, it follows 
that ||re>||c0’i([i,oo)x[o,i]) ^ E. Using Ux = h~^(t)wy again, we get ()2.4p . The proof is complete, q 

3 Preliminaries 

In this section, we first state a comparison principle and then show some properties of the principal 
eigenvalue of a T-periodic eigenvalue problem. Finally, we discuss the existence and uniqueness of 
positive solution to a T-periodic boundary value problem in the bounded interval. 

The following lemma is the analogue of Lemma 3.5 in m and the proof will be omitted. 
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Lemma 3.1 {Comparison principle) Let h G C'^([0,oo)) and /i > 0 in [0,oo), u G C^’^{Q) H 
with Q = {{t,x) : t > 0, 0 < X < h{t)}. Assume that {u,h) satisfies 

ut — duxx > o,{t, x)u — b{t, x)vfi‘, t > 0, 0 < X < h{t), 

< i3[n](t, 0) > 0, u{t,h{t)) = 0, t>0, 

, h'{t) > —fiUx{t, h{t)), t > 0. 

If h{0) > ho, u{0,x) > 0 in [0,h{0)], and u{0,x) > uo{x) in [0,/io]. Then the solution {u,h) of 
dni) satisfies h{t) < h{t) in [0,oo), and u <u in D, where Zl = {(t, x) : t > 0, 0 < x < h{f)}. 

For any given Z > 0, let Xi{£]d,a) be the principal eigenvalue of the T-periodic eigenvalue 
problem 

fit — dfixx — a(L x)fi = I'fi, 0<t<T, 0 < x < i, 

< 0) = 0, fi{t,i) = 0, 0<t<T, (3.1) 

_ fi{0, x) = fi{T, x), 0 < X < 1. 

Proposition 3.1 The principal eigenvalue Ai(Z; d, a) is continuous and strictly decreasing in a and 
£. Moreover, lim£_^o+o) = oo and 

lim Ai(Z; d, a) = oo. (3.2) 

d^oo 

Proof. We only prove (13.21) . the other conclusions can be found in the monograph [20]. Define 
a = supj 2 ,>o o(t, x). Then Ai(f;d, a) > \i{t,d,d) since \i{t,d,a) is decreasing in a. Because a is a 
constant, we know that Xi{£;d,d) is the principal eigenvalue of the elliptic problem 

{ —dfi" — afi = Xfi, 0 < X < I, 

B[fi]{0) = 0, fi{i) = 0, 

and limrf_,.oo Ai(£; d, d) = oo. Thus (|3.2D holds. This completes the proof. q 

In order to study the spreading phenomenon and establish the sharp criteria in later, we shall 
introduce some sets and analyze their properties. For any given d > 0, dehne = {£> 0 : 
Xi(£; d, a) =0}. By the monotonicity of Xi{£] d, a) in £, we see that contains at most one element. 
For any given Z > 0, we let Yle = {d > 0 : Xi{£]d,a) < 0} and Ylt = {d > 0 ■ Xi{£;d,a) > 0}. 
Here we should remark that because Ai(f; d, a) is not monotone in d (cf. Theorem 2.2 of |21]1. it is 
useless to define as the manner treating 
First of all, Ylt 7 ^ ® by Proposition 13.11 

Remark 3.1 For the fixed d > 0, because Xi{£; d, a) = oo and \\m.^^ooXi{£-,d,a) := 

Xi{oo‘,d,a) exists, we have that Yfd 7^ 0 equivalent to Xi{oo]d,a) < 0. 

Proposition 3.2 Assume that the function a{t, x) satisfies 

(HI) There exist ? > 0, —2 < p < Q, k > 1 and Xn satisfying x^ —>■ oo as n ^ oo, such 
that a{t,x) > oxP in [0,T] x [xn,kxn]. 
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Then Ai(oo; d, a) < 0, and so / 0 for any d > 0. 


Proof. Let Xi{kxn',d,a) and 7 f^(xn;d, a) be the principal eigenvalues of 

(ft — d(fxx — CL{t, x)(f = A(/>, 0<t<T, 0<X< kXn, 

(f{t, 0) = 0, (j){t, kxn) =0, 0 < t < T, 


, 0(0, x) = 0(T,x), 


0 < X < kXn 


and 


'ft - dfxx - a{t, x)f = 70, 0 <t <T, Xn < X < kXn, 

0(0 Xn) = 0, 0(t, kXn) =0, 0 <t <T, 


, 0(0, x) = 0(T,x), 


Xn < X < kXn, 


respectively. Then {kxn', d, a) < 'yf^{xn',d,a) by Proposition 13.11 In view of Proposition 17.7 in 
pO] . Xi{kxn',d,a) < Xf{kxn',d,a). Thanks to /? < 0, one has a(t,x) > c^x^ > (^kPxn in [0,T] x 
[xn, kxn]- Since 7f’(xn; d, a) is decreasing in a, we have 7f’(xn; d, a) < 7f’(xn; d, (^k^Xn)- Therefore, 


Xi{kxn]d,a) < yf (xn;d,?Fx^), V n > 1. 


(3.3) 


Let 0(t,x) be the positive eigenfunction corresponding to 7 f’(xn; d, ^/c^Xn). Set y = x/xn and 
T(t, 7/) = 0(t,x). Then ^{t,y) satisfies 

Tt - dxf^'^yy - qkPxfi^ = yf (xn; d, ?/i;^x^)T, 0 <t <T, 1 < y < k, 

< 'I'(01) = O, 'I'(0A;) = O, 0 < t < T, 

_ T(0,y) = T(r,7/), l<y<k. 

Utilizing the inequality (17.6) of [20], we have 

7f’(xn; d, (ik^xf) = dX*x~‘^ — ?Kx^ = x“^(dA* — (ik^x^P) <0 for 1 (3.4) 

since 2 + p > ^ and Xn —)• oo as n ^ oo, where A* is the principal eigenvalue of 

—u” = Xu, 1 < y <k, 
u(l) = 0, u{h) = 0. 

It follows from (13.3p and (13.4p that Ai(oo;d, a) < 0. The proof is complete. q 

The condition (HI) seems to be “weak” because a{t, x) may be “very negative” in the sense 
that both |{a(t,x) > 0}| <C |{a(t, x) < 0}| and a(t, x)dx = —oo are allowed for any t > 0. 


Proposition 3.3 Assume that 

(H2) there exists x > 0 such that a(t, x}dt > 0. 
Then for any £ > x, the set is non-empty. 
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Proof. Owing to a{t, x)dt > 0, there exists 0 < e < min{x, £ — x} such that a{t)dt > 0, 
where a{t) = min^jg/ a{t, x) and I := [x — e,x + e\. Let A be the principal eigenvalue of 

{ —u” = \u, X — e<x<x + e, 
u{x ± e) = 0. 

Then A > 0. Applying the inequality (17.6) of [5D] we have 

Xi{£', d, a) < dX — — / a{t)dt. 

Jo 

rp 

Owing to Jq gJ(t)dt > 0, there exists do > 0 such that Xi(£;d,a) < 0 for all 0 < d < do- This 

implies ( 0 , do) C □ 

Now we consider the following T-periodic boundary value problem of logistic equation in a 
bounded interval (0,f): 

vt — dvxx = a(L x)v — b{t, x)v‘^, 0 < t <T, 0 < x < £, 

< 0) = 0, v{t,£) = 9, 0<t<T, (3.5) 

u(0, x) = v{T, x), 0 < X < £. 


Lemma 3.2 Assume that the function a satisfies (HI) and Iq is the unique positive root of 
Xi{i]d,a) = 0. Then for any given £ > £o and 9 > ||a||oo/^i, the problem (13.5p has a unique 
positive solution, where bi is given by (H). 


Proof. The approach used in this proof can be regarded as the upper and lower solutions 
method. Since £ > £o, we have Xi{£]d,a) < 0. Let (j) be the positive eigenfunction of (13.111 
corresponding to Ai(f;d, a). It is easy to verify that ecj) is a lower solution of (13.511 and efi < 9 
provided 0 < e ^ 1. 

Let z be the unique solution of the initial-boundary value problem 

zt — dzxx = a(L x)z — b{t, x)z‘^, t > 0, Q < x < £, 

< B[z]{t,0) = 0, z{t,£) = 9, t > 0, 

z(0, x) = 9, 0 < X < £. 

Then £((){£, x) < z{t, x) < 9 by the comparison principle for parabolic equations. For the integer 

n > 0, we define z^{t,x) = z{t + nT,x), {t,x) G [0,T] x [0,.^]. Because a,b are T-periodic in t, we 
see that z” satisfies 

zf' — dz^x; — ®(^) x)z^ — b{t, x)(z^)‘^, 0 < t <T, 0 < x < £, 

< B[z^]{t,0) = 0, z^{t,£) = 9, 0<t<T, 

_ z'^(0, x) = z(nT, x), 0 < X < £. 

Since 

£c/){0,x) = e(f{T,x) < z{T,x) = z^(0,x), z^(0,x) = z{T,x) <9 = z(0,x), 
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we can apply the comparison principle to get that ecf) < < z in [0,T] x [0,£]. Which also implies 

60(0,x) < z^{i),x) = z^{T,x) < z{T,x) = z^(0,x). 


As above, ecj) < z^ < z^ in [0,T] x [0,£]. Applying the inductive method we have that z^ is 
decreasing in n and z^ > scj) in [0,T] x [0,£]. So, there exists a non-negative function v such 
that z"’ —^ u pointwise in [0,T] x [0,i] as n ^ oo. Since z”'+^(0,x) = z"‘{T,x), it follows that 
u(0, x) = v{T,x). Based on the regularity theory for parabolic equations and compact argument, 
it can be proved that there exists a subsequence {ui}, such that z"^* —>■ U in C'^’^([0,T] x [0,^]) as 
i oo, X > 60 in [0, r] x [0,£] and v satisfies the first two equations of (I3.5p . This shows that v is 
a positive solution of (13.Sp . 

Now we prove that v is the unique positive solution of (I3.5I1 . Let v be another one. By the 
maximum principle, we have v{t,x) < 6, thereby, v < z. This implies x < z” in [0,T] x [0,£] for any 
integer n > 0. Certainly, x < x in [0,T] x [0,£]. This suggests that x is a maximal positive solution 
of (13.5p . It is easy to see that there exists a constant 6 > 0 such that x > 6x in [0,T] x [0,^], and 
hence the infimum 


a = 


int 060 

0<t<T v(t,x) 

o<x<e ^ ' 


exists and is positive. Clearly, a < 1 and x > av in [0,T] x [0,£]. If we can show <7 = 1 then x = x 
and the uniqueness is derived. Assume on the contrary that u < 1. Denote (p = v — av, then </? > 0, 
and 


(p{0, x) = <p{T, x) for 0 < X < £, 

B[{p]{t,0) = 0, = {1 — a)6 > 0 for 0<t<T. 

The direct calculation yields 


Pt - dpxx = «(L x)ip - b{t, x) (x^ - cjx^) 

> a{t, x)ip — b{t, x) (x^ — cj^x^) 

> a{t, x)ip — 2b{t, x)vip. 

The maximum principle allows us to deduce that (^ > 0 in [0,T] x (0,£]. Similarly to the above, 
there exists 6i > 0 such that p > 6ix, and thus x > {a + ei)v in [0, T] x [0, i]. This is a contradiction 
with the definition of a, and the uniqueness is derived. q 

4 Spreading-vanishing dichotomy and long time behavior of solu¬ 
tion 

We first give a lemma, by which the vanishing phenomenon is immediately obtained. Moreover, 
this lemma will play an important role in the establishment of sharp criteria for spreading and 
vanishing. 
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Lemma 4.1 Let d, g, and B be as above, c £ M.. Assume that g £ C^{[0,oo)), ip £ C 2 00 ) x 


[0,5((i)]) and satisfy g{t) > 0, ip(t,x) > 0 for t > 0 and 0 < x < g{t). We further suppose that 
limt_>.oo 5(t) < 00 , = 0 and there exists a constant C > 0 such that •)llci[o,g(t)] ^ 

C fort > 1. If [ip, g) satisfies 


/ 


t > 0, 0 < X < g{t), 

i > 0, X = 0, 

t>0, x = g{t), 


P't - dipxx > op, 

< B[ip] = 0, 

^ = 0, g'{t) > -gipx 


then limt_^oo maxo<j,<g(t) ip{t, x) = 0. 

Proof. When a = 0 or /3 = 0, this is exactly Proposition 3.1 of [35]. When a > 0 and /3 > 0, 
that proof remains valid. The details are omitted here. q 

If /loo < 00 , the estimate (12.3p implies limt_).oo h'{t) = 0. Applying (|2.3p . ()2.4I1 and Lemma ITTl 
we have the following theorem. 

Theorem 4.1 {Vanishing) Let {u,h) be the solution of (11.11) . When h^o < 00 , we must have 


lim max u{t, x) = 0. 

t^oo 0<x<h(t) 


(4.1) 


This shows that if the species cannot spread successfully, it will be extinct in the long run. 


In the following we investigate the spreading phenomenon. To this aim, we first study the 
existence and uniqueness of positive solution to the following T-periodic problem 


Ut — dUxx = o.{t, x)U — b{t, x)U‘^, 0 < t < T, 0 < x < 00 , 


< B[U]{t,0)=0, 
^Ui0,x) = U{T,x) 


0<t<T, 

0 < X < 00 . 


(4.2) 


Theorem 4.2 Assume that there exist a constant p with —2 < p < 0, and T-periodic positive 
functions aoo{t),boo{t),a°°{t),b°°{t) £ C^/^([0, T]), such that 


/ 


dooit) = lim inf , a°°{t) = lim sup 

x^co xP X^OO 



(4.3) 


boo{t) = liminf 6(/, x), 6°°(/) = limsup6(/,x) 



uniformly on [0,T]. Then (|4.2I) has a unique positive solution U £ C'^''“2’^+'^([0, T] x [0,oo)), and 
satisfies 



(4.4) 


uniformly on [0,r]. 


It is worth mentioning that the condition (031) implies the assumptions (HI) and (H2). 
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Proof. When a = 0, this theorem is a special case of Theorem 1.3 in [32]. We consider the 
case a > 0 in the following. This proof is divided into three steps. In the first one, we construct 
the minimal positive solution of (14.21) . The estimate ()4.4p will be given in the second step. Finally, 
we show the uniqueness of positive solution. 

Step 1: The existence. In this step, we shall construct a positive solution U_ and prove that it 
is the minimal one. Consider the following auxiliary problem 

Ut — dUxx = «(T x)U — b{t, x)U‘^, 0<t<T, 0<x<f, 

< B[17](t,0) = = 0, 0<t<T, (4.5) 

^ 17(0, x) = U{T,x), 0 <x < i. 

Let Xi{i]d,a) be the principal eigenvalue of (13.ip . Since the assumption (HI) holds, by use of 
Proposition 13.21 there exists £o ^ 1 such that Xi{i,d,a) < 0 for all i > Iq. For such i, utilizing 
Theorem 28.1 of [20], the problem (14.51) admits a unique positive solution, denoted by Ui{t,x). 
Obviously, Ug, < ||a||oo/^i by the maximum principle, where bi is given by (H). For I* > it is 
easy to see that Ug* is an upper solution of (|4.5p . Let 4‘{t,x) be the positive eigenfunction of (13.11) 
corresponding to Xi{i,d,a) and e > 0 be a constant. Then ecj) is a positive lower solution of (14.51) 
and ecj) < Ug* provided e <C 1. Thus, Ug* > Ug since Ug is the unique positive solution of (14.51) . 
This shows that Ug is increasing in Make use of the regularity theory for parabolic equations 
and compact argument, it can be proved that there exists a subsequence of {Ug}., denoted by itself, 
and a positive function U_ G C'^’^([0,T] x [0,oo)), such that Ug ^ U_m. C'^’^([0,T] x [0,Lj) for any 
L > 0, and U_ solves 

Let 17 be a positive solution of (|4.2I) . Then U < ||a||oo/l'i by the maximum principle. Obviously, 
U is an upper solution of (14.51) for any given £ > Iq. As above, Ug < U in [0, T] x [0 ,£]. This implies 
U_<U< ||a||oo/fti in [0,7"] x [0,oo). 

Step 2: Proof of ()4.4I) . For any positive solution U of ()4.2p . we have known that U_ < U < 
||o||cxd/1)i in [0,T] X [0, oo). Take account of Theorem 1.3 in [52], the problem 

Wt — dWxx = o(L x)W — b{t, x)W‘^, 0 < t <T, 0 < x < oo, 

< Wxit,0) = 0, 0<t<T, (4.6) 

_ W{0,x) = W{T,x), 0 < X < oo 

has a unique positive solution W{t,x). For £ > 0, let us consider the problem 

Vt — dVxx = n(L x)V — b{t, x)I7^, 0 <t <T, 0 < x < £, 

< Vx(t,0) = 0, V{t,£) = U{t,£), 0<t<T, (4.7) 

^V{0,x) = V{T,x), 0<x<£. 

It is obvious that U is a lower solution of (|4.7p since Ux{t, 0) > 0, and the constant K = 1 + ||o||oo/l)i 
is an upper solution of (|4.7I) . Similarly to the proof of Lemma 13.21 we can prove that the problem 
(iZI) has a unique positive solution Vg and U < Vg < K in [0,T] x [0,£]. 

Arguing as Step 1, there exists a subsequence of {V^}, denoted by itself, such that Vg ^ W 
in C'^’^([0,T] X [0, L]) for any L > 0 since W is the unique positive solution of (14.61) . Therefore, 
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U <W in [0,T] X [0, oo). In view of Eq. (2.6) in [32], we have that 

U{t,x) W{t,x) Tnax[o,T]a°°{t) 

hm sup - < hm sup - < -;- 

x^oo xP x^oo xP min[0,r] Ooo(i) 


uniformly on [0,r]. 
Set 


a(x) = mina(t, x), 6(x) = max6(t, x), 
[o,T] [o,r] 


and consider the problem 


—dw” = a{x)w — b{x)w‘^, 0 < x < oo, 

t(;(0) = 0. 


(4.8) 


(4.9) 


Thanks to the conditions (H) and g3D, we can show that a(x), b{x) G (^^([0, oo)), and 


0 < minaoo(i) = liminf < limsup < max. a°°(t), 
[0,T] x^oo xP x^oo XP [0,r] 

0 < min6oo(t) < liminf 6(x) < limsup6(x) = max6°°(t). 
[ 0 ,T] x^oo x^oo [ 0 ,T] 

By virtue of Proposition 13.21 the first eigenvalue of 


— a{x)(f> = \4>, 0 < x < £, 

m = o, m = o 


is negative provided i ^ 1. And then, the problem 


—dw" = a{x)w — b{x)w‘^, 0 < x < i, 

u;(0) = 0, wii) = 0 


has a unique positive solution Wi{x). By the same argument as in Step 1, we can show that there 
exists a positive function w such that wi ^ w in C'^’^([0,L]) for any L > 0, and w solves (|4.9p . 
Obviously, wi is a lower solution of (|4.5I) . Hence, wi < since is the unique positive solution 
of (14.5p . This implies w < U_ < U in [0,T] x [0, oo) since wi ^ w and Ui —>■ IJ_. Moreover, t(;'(0) > 0 
by the uniqueness. 

If w'{x) > 0 in [0, oo), then w{x) ^ tc* as x —>■ oo for some positive constant tc*. As /? < 0, it 
is immediate to get 


x^oo x^oo 


uniformly on [0, T]. In this case we claim that p = 0 and 

^ lim inf3;^oo a(a^) ™™[o,T] ®oo(i) 

limsup,,^^6(x) max[o,r] b°°{t) 


(4.10) 


(4.11) 


In fact, if p < 0, then lima,^oo a(^) 3^) = 0 uniformly on [0,r] by (14.31) . Therefore, lim 3 ;_j.oo 0 ( 2 ;) = 0. 
This is impossible since b{x) has a positive lower bound, w'{x) > 0 in [0,oo) and w is bounded 
from above. So, p = 0. On the contrary we assume that (14.111) does not hold. Denote a* = 
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lim inf 3 ;_>.oo a(x) and b* = limsupa,_^oo b{x). Then there exist e > 0 and x* S> 1, such that w{x) < 
a^f/b* — £ for all x > x*. For such e > 0, there exist 5 > 0 and xq > x* such that 

- S a* 

b* + 6 ~ b* ^ ^ ^ ^ ’ 

and a(x) > — 5, b{x) < b* + 6 for all x > xq. Then 

—dw" = w[a{x) — b{x)w] > wb{x) f ^ ^ ) > crb{x)w, V x > xq, 

V Kx) 0* / 

this is impossible since cr > 0 and 5(x) has a positive lower bound. So, (14.lip is true. 

If w'{xo) = 0 for some xq > 0, we set 

a*(x) = a(x + Xq), 6*(x) = 6(x + Xq), u(x) = u;(x-|-Xq) for X > 0, 

a*(x) = a(— X -|-Xq), 6*(x) = 6(— X -|-Xq), v{x) = w{—X + Xq) for X < 0. 

Then the function v{x) satisfies 

— dv" = g*{x)v — b*{x)v‘^, —oo < x < oo. (4.12) 

Applying Theorem 7.12 of [8], we have that v{x) is the unique positive solution of (|4.12p and 
satisfies 

a(x) 


liniinf^> ™ 


lim inf ^ 

X^OO 


mm 


[0,T] Ooo(i) 


Therefore, 


xP lim sup 6* (x) limsup6(x) max[o,T] ' 

a:—>-co x^oo 


lim inf Idlh > lim inf = lim inf ““<*> 


XP 


(4.13) 


x^oo xP x^ao xP maXjo^T] b°°{t) 

uniformly on [0,r]. 

It follows from (14.81) . (|4.10p . ()4.1ip and ()4.13p that U satisfies (|4.4I) . 

Step 3: The uniqueness. Let U he a positive solution of (|4.2I) . Then U > U_, here IJ_ is the 
minimal positive solution of (14.211 obtained by Step 1. To prove the uniqueness, it suffices to show 
that U = IJ_. If this is not true, then U ^ U_. It follows from (j4.4l) that there exists k > 1 such 
that U < kU_ in [0,T] x [0,oo). To arrive at a contradiction, we turn to a technique introduced by 
Marcus and Veron in [30]. Define V = U_— {2k)~^{U — U_). Then 


c, > > itij;, 


2k 


-fo + 


1 


-U = U. 


(4.14) 


2k — 2k+ 1 2k+ 1 

Noticing that is convex in z G (0,oo), we have second formula of 

(|4.14p . Because b{t,x) is positive, then the direct computation gives 

Vt-dV,,>a{t,x)V-b{t,x)V\ 

Obviously, for the large £, 

B[V]{t,0) = 0, V{t,i)>0 in [0,r]; V{0,x) =V{T,x) in [0,£]. 
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This shows that V is an upper solution of (I4.5p . Let [/£ be the unique positive solution of (14.51) . As 
above, U'£ < V in [0,T] x [0,.^]. Hence, ^ H in [0, T] x [0, oo) due to ^ as £ —>■ oo. This is 
a contradiction with the first inequality of (I4.14D . So, U = U_ and the uniqueness is derived. The 
proof is complete. q 

Theorem 4.3 (Spreading) Assume that (j4.3l) holds. If h^o = oo, then lim^^oo 'U'(t + x) = 
U(t,x) uniformly in [0, T] x [0,L] for any L > 0. 

Proof. At present, the condition (HI) holds. Let Xi(i;d,a) be the principal eigenvalue of 
(lOl . In view of Proposition 13.21 we can choose an fo 1 such that Xi{£,d,a) < 0 for all i > io- 
We divide the proof into two parts, and show that, respectively, 

limsup u{t + nT,x) < U(t,x) uniformly in [0,T]x[0,L] (4.15) 

n—^oo 

and 

liminiu{t + nT,x) > U{t,x) uniformly in [0,T]x[0,L]. (4-16) 

71—^00 

Step 1. Take 9 = ||no||oo + ^||a||oo 5 where hi is given by (H). Then u < 9 hy the maximum 
principle. For the fixed I > (.q, there exists an integer m 3> 1 such that h(t) > I for all t > mT. 

Since Xi{i]d, ||a||cxD) < Xi{i-,d,a) < 0, by the upper and lower solutions method we can prove 
that the boundary value problem 

j-dVxx = \\a\\ooV-biV‘^, 0 < x < £, 

|h[P](0) = 0, V(i) = 9 

has a unique positive solution V(x). In consideration of the regularity of u(mT,x) and V{x) in 
X, we can find a constant k > 1 such that u{mT,x) < kV(x) for all 0 < x < f. Note that 
u(t,i) < k9 = kV(i), we can apply the comparison principle to u and kV, and then derive that 
u(t,x) < kV(x) for all t > mT and < x < L Since /c > 1, it is easy to see that the function 
V := kV satisfies 

—dvxx > llolloo'i^ — biv"^ > a{t, x)v — b(t, x)v^. 

Let W£ be the unique solution of 

wt — dwxx = o,(t, x)w — b{t, x)uP‘, t > mT, 0 < x < f, 

< H[r(;](t, 0) = 0, w{t,i) = kV(i), t > mT, 

_ w(mT, x) = kV (x), 0 < X < £. 

The comparison principle gives u <W£ < kV for t > mT and 0 < x < f. For the integer n > m, we 
define Wg{t,x) = W£(t + nT,x), (t,x) G [0,r] x [0,f]. Similarly to the proof of Lemma it can 
be shown that W£ in C^’^([0,r] x [0,f]) as n ^ oo, where W£ is the unique positive solution 

of (j3.5|) with 9 replaced by k9. Owing to u(t + nT,x) < W£(t + nT,x) = w^{t,x) in [0,r] x [0,.^], 
we get 

limsupn(t + nP, x) < W£{t,x) uniformly on [0, T] x [0,f]. 

n—>-co 


(4.17) 
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Since < k6, it is easy to see that is decreasing in i. Remember that U is the unique 
positive solution of (14.21) . arguing as Step 1 in the proof of Theorem 14.21 we can prove that U 

in X [0,L]) as £ — oo. Combining this fact with (j4.17p . we get ()4.15p . 

Step 2. Let £ > io and vi be the unique solution of 

vt — dvxx = 0'{t, x)v — b{t, x)v‘^, t > mT, 0 < x < £, 

< 0) = 0, v{t,£) = 0, t > mT, (4.18) 

_ v{mT,x) = u{mT,x), x € [0,£]. 

Then u > ve in [mT,oo) x [0,^]. Since \i{£;d,a) < 0, the problem 

Vt — dVxx = a{t, x)V — b{t, x)V'^, 0 < t <T, 0 < x < £, 

< 0) = 0, V{t,£) = 0, 0<t<T, 

^V{0,x) = V{T,x), 0 < x < £ 

admits a uniqe positive solution V^. Clearly, the constant 6 = ||uo||oo + ^ ||fl||oo is an upper solution 
of (14.181) . Noticing Vi{mT,x) = u{mT,x) > 0 in (0,£], take advantage of Theorem 28.1 in [20], it 
follows that Vi{t + nT,x) —>■ Vi(t,x) in C'^’^([0,T] x [0,£]) as n —>■ oo. 

It is obvious that Vi is increasing in £. Similarly to Step 1, we can derive that Vi ^ U in 
X [0, Lj). Consequently, the limit (|4.16l) is obtained since u > vi in [mT, oo) x [0,£] for 
all £ > £o. The proof is finished. □ 

5 Sharp criteria for spreading and vanishing, spreading speed 

We first give a necessary condition for vanishing. Let Xi{£', d, a) be the principal eigenvalue of (13.ip . 

Lemma 5.1 If hoo < oo, then Xi{hoo',d,a) > 0. 

Proof. We assume Xi{hoo',d, a) < 0 to get a contradiction. By the continuity of Ai(£; d, a) in £ 
and h{t) —)• hoo, there exists r S> 1 such that Ai(/i(r); d, a) < 0. Let w be the unique solution of 

Wt — dwxx = o(t, x)w — b{t, x)wP‘, t > T, 0 < X < h{T), 

< B[w]{t,0) = w{t,h{T)) = 0, t > T, 

_ w{t,x) = u{t,x), 0 < X < /i(r). 

Then n > re in [r, oo) x [0 ,/i(t)] by the comparison principle. Remembering Ai (/i(t) ; d, a) < 0, 
it follows from Theorem 28.1 of [20| that w[t + nT,x) z{t,x) uniformly in [0,T] x [0 ,/i(t)] as 
n —>■ oo, where z is the unique positive solution of the following T-periodic boundary value problem 

zt — dzxx = o,{t, x)z — b{t, x)z‘^, 0 < t <T, 0 < x < h{T), 

< B[z]{t,0) = z{t,h{T)) =0, 0 < t < T, 

_ z(0,x) = z{T,x), 0 < X < /i(r). 

Since u > w in [r, oo) x [0,/i(r)], it is deduced immediately that 

liminf u{t + nT, x) > z{t, x), V (t, x) £ [0, T] x [0, ^(t)]. 

n—>-oo 

This is a contradiction with (14.ip . The proof is complete. q 
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Lemma 5.2 If Xi{ho;d, a) > 0, then there exists go > 0, such that h^o < oo provided g < go- 
Hence, by Lemma \5A\ \i{hoo',d,a) > 0 for g< go- 


Proof. Let 4>{t, x) be the corresponding positive eigenfunction to Ai := Ai(/io; d, a) of (13.11) with 
£ = ho- Noticing that (fxit, ho) < 0, (p{t, 0) > 0 in [0, T] when (3 > 0, and 0) > 0 in [0, T] when 

/3 = 0. By the regularity of (f, there exists a constant C > 0 such that 


x(f>x{t,x) < C4>{t,x), V {t,x) G [0,T] X [0, ho]. (5.1) 

Let 0 < (5, (T < 1 and K > 0 he constants, which will be determined later. Set 

s{t) = 1 + 2(5 — 6e~^*, T{t) = f s~‘^{p)dp, t > 0, 

Jo 

v{t,x) = Ke~‘^^(l){T{t), y), y = y{t,x) = -^, 0<x<hos{t). 

Firstly, for any given 0 < e <C 1, since a is uniformly continuous in [0, T] x [0, S/iq] and T-periodic 
in t, we have that there exists 0 < 5o(e) ^ 1 such that, for all 0 < J < 6o{e) and 0 < cr < 1, 


|s ‘^{t)a{T{t), y{t,x)) — a{t,x)\ <e, Vt>0, 0<a:< hos{t). 

Note that (15.Ih . (15.2j) and Ai > 0, the direct calculation yields, 

J N f , a{T,y) (. ^ y4>y{r,y) crd Ai 

vt — dvxx — a.{t, x)v = V I —a H—- a,{t, x) -77-7- —e + 


(5.2) 


s2(t) ’ </>(r,y) s{t) s‘^{t) 

> v{—(j — £ — Ca + Ai/4) >0, Vt>0, 0<x< hos{t) 


(5.3) 


provided 0 < cr, e <C 1. 

Evidently, v{t,hos{t)) = Ke~^^f){T{t),ho) = 0. If either a = 0 or /3 = 0, then i?[u](t,0) = 0. If 
a,/3 > 0, then a(p{T{t),0) = j3(f>y{T{t),0), and (j)y{T{t),0) > 0 by the Hopf lemma. Therefore, 

5[u](t,0) = f3Ke~'"^4>y{T{t),0)[l - l/s(t)] > 0 

due to s(t) > 1. In a word, 

i?[u](t, 0) > 0, v{t,hos{t)) = 0, y t > 0- (5.4) 

Fix 0 < CT, e <C 1 and 0 < <5 < 6o{£)- Based on the regularities of uo{x) and (/>(x), we can choose a 
constant K ^ 1 such that 


^o(a^) <(0, x/(l + 5)) = u(0, x), y 0 < X < ho- (5.5) 

Thanks to hos'{t) = hocr6e~'^^ and Vx{t,hos{t)) = ^iFe“°'Vy(r(t), ho), there exists go > 0 such 
that, for all 0 < ^ < go, 

hos'{t) >—gVx{t,hos{t)), V t > 0. (5-6) 

Remember (I5.3p - (|5.6p . Applying Lemma [3T] to {u,h(t)) and {v,hos{t)), it follows that 
h{t) < hos{t), u{t,x) < v{t,x), Vt>0, 0<x< h(t). 

Hence h^o < hos(oo) = ho(l + 26) for all 0 < /U < go- This finishes the proof. q 
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Lemma 5.3 Let C > 0 be a constant. For any given constants ho,H > 0, and any function 
uq E /iq]) satisfying 71[uo](0) = uo(/io) = 0 uq > t) in (0,/io); there exists /x° > 0 such 

that when fj,'> fjP and {u, h) satisfies 


Ut — duxx > —Cu, 

B\u]{t, 0) = 0 = u{t, h{t)), 
h'{t) = -^Ux{t,h{t)), 


t > 0, 0 < a: < h{t), 
t > 0, 
t > 0, 


_ u{0,x) = uq { x ), h{0) = ho, 0 < X < /iQ, 


we must have liminft^oo h{t) > H. 


The proof of Lemma (5.31 is essentially similar to that of Lemma 3.2 in [38] and is hence omitted. 
Now we fix d, and consider ho and ^ as varying parameters to depict the sharp criteria for 
spreading and vanishing. Assume that / 0 and let h* = h*{d) E ’Yhd-' ^i{h*',d,a) = 0. It 
is worth mentioning that if the assumption (HI) holds, then / 0. 

Recalling the estimate ()2.2p , as the consequence of Lemmas 15.1115.21 and 15.31 we have 


Corollary 5.1 (i) If hoo < oo, then hoo < h*. Hence, ho > h* implies hoo = oo for all /r > 0; 

(ii) When ho < h*, there exist pLo, > 0, such that hoo < h* for p, < po> ond hoc = oo for 

Finally, we give the sharp criteria for spreading and vanishing. 


Theorem 5.1 (i) If ho > h*, then hoc = oo for all p > 0; 

(ii) If ho < h*, then there exists p* > 0, such that hoc = oo for p > p*, while hoc < h* for 
p<p*. 

Proof. Noticing Corollary 15.11 by use of Lemma l3.1l and the continuity method, we can prove 
Theorem O Please refer to the proof of Theorem 3.9 in m for details. □ 

Now we fix ho, and regard d and p as the variable parameters to describe the sharp criteria for 
spreading and vanishing. First of all, Ylto 7^ ® Proposition 13.11 

Theorem 5.2 (i) When d E have hoc = oo for all p > 0; 

(ii) For any fixed d E Ylho’ I'^ore exists po = Po{d) > 0 such that hoc < oo provided < p < po- 
If, in addition, Yld 7^ ® Z®’" such d, then there exists p* > 0, such that hoc = oo when p > p*, and 
hoc < oo when p < p*. 

Remark 5.1 (i) Let x be given in (H2). From the proof of Proposition 13.31 we see that if ho > x 
then 7 ^ 0, and there exists do > 0 such that d E all 0 < d < do- 

(ii) If the condition (HI) holds, then Jf^d 7^ ® Z®’" all d > 0 by Proposition 13.21 

Proof of Theorem 15.21 (i) For any d E have \i{ho]d,a) < 0. If \i{ho',d,a) < 0, 

then Yld 7^ ® ^0 > h*{d). If Ai(/io; d, a) = 0, then ho = h*{d). By Theorem l5.11 if . hoc = oo for 

all p > 0. 
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(ii) For the fixed d G have Ai(/io;d, a) > 0. By Lemma [5.21 there exists go > 0 

such that /loo < oo for /X < gQ. If, in addition, Yld then there exists iL S> 1 such that 

Xi{H-,d,a) < 0. Take advantage of Lemma [531 there exists g^ > 0 such that hoo > H provided 
g > g^i which implies Xi{hoo]d,a) < Xi{H;d,a) < 0. Hence, hoo = oo for /x > /x^ by Lemma [5Tl 
The remaining proof is the same as that of Theorem 3.9 in m- □ 

Now we fix d and consider ho,g as the varying parameters. Note that (14.3p implies (HI), 
combining Theorems 14.11 14.31 and 15.11 we immediately obtain the following spreading-vanishing 
dichotomy and sharp criteria for spreading and vanishing. 

Theorem 5.3 Assume that (j4.3p holds. Let h* be the unique positive root of Xi{i;d, a) = 0 and 
{u,h) be the unique solution of (II.ip . Then the following alternative holds: 

Either 

(i) Spreading: hoo = oo and lim„_j.oo «(/ +nT, x) = U{t,x) uniformly on [0,T] x [0, L] for any 
L > 0, where U{t,x) is the unique T-periodic positive solution to (14.21) .- 

or 

(ii) Vanishing: hoo < h* and lim^^oo ^^'^o<x<h{t) u{t, x) = 0. 

Moreover, 

(iii) If ho > h*, then hoo = oo for all g > 0; 

(iv) If ho < h*, then there exist g* > 0, such that hoo = oo for g > g*, while hoo < h* for 
g<g*. 

When we fix ho and regard d, g as the variable parameters. Note that (|4.3p implies (HI), using 
Theorems [mono and Remark 15. 11 we have the following spreading-vanishing dichotomy and 
sharp criteria. 

Theorem 5.4 Assume that ho > x, where x is given in (H2). Let the condition (|4.3p hold and 
{u,h) be the unique solution of (jl.ll) . Then the following alternative holds: 

Either 

(i) Spreading: hoo = oo and lim„_).oo n(t-|-riT, x) = U{t,x) uniformly on [0,T] x [0, L] for any 
L > 0, where U{t,x) is the unique T-periodic positive solution to (j4.2p .- 

or 

(ii) Vanishing: hoo < oo and limi_).oo maxo<a,</i(t) u{t,x) = 0. 

Moreover, 

(iii) If d G Ylho’ ~ h ^ 

(iv) If d G ^ho’ then there exists g* > 0, such that hoo = oo when g > g*, and hoo < oo when 
g<g*. 

In the last part of this section, we estimate the asymptotic spreading speed of the free boundary 
h{t) when spreading occurs. To this aim, let us first state a known result, which plays an important 
role in the study of asymptotic spreading speed. For an integrable T-periodic function f{t), we define 

- 1 r 
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Proposition 5.1 f [111 Section 2]) Let d > 0 and 0 < < 1 be the given constants. Assume that 

p,q £ C'‘^([0,T]) are positive T-periodic functions and k G C'‘^([0,T]) is a nonnegative T-periodic 
function. Then the problem 

wt — dwxx + k{t)wx = p{t)w — q{t)uP‘, 0 < t < T, 0 < x < oo, 

< w{t, 0) = 0, 0 < t < T, 

^ w{0, x) = w{T, x), 0 < X < oo 

has a positive T-periodic solution w^{t,x) € C^'’^{[Q,T] x [0,oo)) if and only if k < 2y/^, and such 
a solution is unique when it exists. Furthermore, the following hold: 

(i) w^{t,x) > 0 and w^{t,x) v{t) uniformly on [0,T] as x ^ oo, where v{t) is the unique 
positive periodic solution of the problem 

v'= p{t)v — q{t)v‘^, 0<t<T; v{0) = v{T); 

(ii) For any given nonnegative T-periodic function m G C'^{[0,T]) satisfying m < 2^/dp, the 
assumption m <, ^ k implies wff{t,0) > Wx{t,0), w'^{t,x) > w^{t,x) for 0 <t <T and x > 0; 

(iii) For each p > 0, there exists a unique positive T-periodic function ko{t) = ko{p,p,q){t) G 
C^dOjT]) such that pWx°{t,0) = ko{t) in [0, T], and 0 < ko < 2\/dp. 

Theorem 5.5 Assume that (14.3p holds with p = 0. When the spreading occurs, i.e., hoo = oo, we 
have 

liminf / koip,aoo,b°°){t)dt, limsup / koip,a°°,boo){t)dt, 

t I Jo t^oc t I Jo 

Proof. This proof can be done by the same manner as section 4 of P. Because the length is 
too long, we omit the details. The interested readers can refer to that reference. 

6 Concluding remarks 

We comment finally on some points raised by the theoretical investigation. Firstly, from the above 
discussion we have seen that Ai(oo;d, a) := lim£_).oo Ai(^; d, a) < 0 is an essential condition. This 
number is only characterized by the dispersal rate d and intrinsic growth rate a{t, x), is independent 
of the self-limitation coefficient b{t,x), moving parameter p and initial value uo{x). 

The main conclusions of this paper can be briefly summarized as follows: 

(I) If one of the following holds: 

(i) the diffusion rate d G Ylho (^Oi “ are fixed); 

(ii) the intrinsic growth rate a{t,x) is suitable “positive” in the sense of Xi{ho]d,a) < 0 {ho, d 
are fixed); 

(iii) the initial habitat ho is suitable “larger” in the sense of Xi{ho',d,a) < 0 {d, a are fixed), 
then the new or invasive species will successfully spread and survive in the new environment (main¬ 
tain a positive density distribution), regardless of initial population size uo{x) and value of the 
moving parameter p. 





20 


Mingxin Wang 


(II) When the above situations do not appear, we can find a critical g,* such that the species 
will spread successfully when g > g*, and that the species fails to establish and will be extinct 
in a long run when g < g*. A better way to modulate the moving parameter g is to control the 
surrounding environment. 

Now we analyze the relations and differences of the free boundary problem dLH) and the corre¬ 
sponding Cauchy problem and fixed domain problem, and make a comparison between our results 
and those obtained for the corresponding Cauchy problem and initial-boundary value problem with 
fixed boundary. 

By the same arguments as that of Theorem 5.4 in m, we can prove that the solution u{t, x) 
of (jl.ip converges to v{t,x) as /r —)• oo, where v{t,x) is the solution of the initial-boundary value 
problem in the half space 


vt — dvxx = a{t, x)v — b{t, x)v‘^, t > 0, x > 0, 
B[v]{t, 0) = 0, t > 0, 

^ v{0,x) = vo{x), x>0, 


( 6 . 1 ) 


where t'o(x) = uo{x) for 0 < x < /iq and vo{x) = 0 for x > /iq. This shows that (16.11) is the limiting 
problem of (jl.ip as /i ^ oo. If a = 0, the problem (16.11) can be extended into the Cauchy problem 


Vt — dvxx = a(T x)v — b{t, x)v‘^, t > 0, x E 
x(0, x) = xo(x), X E M. 


( 6 . 2 ) 


Solutions of both (j6.ip and (16.2p are positive for all x once t is positive. This seems to be a defect 
because the movement should be finite for any species. 

When a and b are positive constants, the differential equation in (16.2p is 


ut — duxx = 'u(a — bu), t > 0, x E 


(6.3) 


A great deal of previous mathematical investigation on the spreading of population has been based 
on the traveling wave fronts and asymptotic spreading speed of ()6.3p . please refer to, for example 
Aronson & Weinberger mm for the model (16.3p . Lewis et al. [26| and Liang &: Zhao [27| for the 
more general models. The known result for (|6.3I) predicts successful spreading and establishment 
of the new species with any nontrivial initial population m( 0, x), regardless of its initial size and 
supporting area. However, this is not supported by empirical evidences, which suggest, in the 
contrary, that success of spreading depends on the initial size of the population; for example, the 
introduction of several bird species from Europe to North America in the 1900s was successful only 
after many or several initial attempts (cf. |34l [29] , where more examples can be found). 

On the other hand, when g = 0, our free boundary problem (II.ip reduces to the following 
initial-boundary value problem with fixed boundary: 

wt — dwxx = a(L x)w — b{t, x)w‘^, t > 0, 0 < x < ho, 

B[w]{t, 0) = w{t, ho) = 0, t > 0, 


^ w{0,x) = uo{x), 


0 < X < ho 


(6.4) 
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These discussions indicate that the free boundary problem dni) is in the intermediate state of (16.111 
and (|6.4p . 

It is well known that (see Theorems 24.2 and 28.1 of [20]) if Ai(/io;d, a) > 0, then the solution 
w of (16.4p must extinct in the long run, i.e., w{t,x) tends to zero as t —)• oo no matter how large 
the initial data uo(x) is and how small the self-limitation coefficient b(t,x) is. This seems do not 
tally with the natural phenomenon. 

For our free boundary problem O), the conclusions stated in Theorems 15.31 and 15.41 show that 

(i) if Xi{hQ]d,a) < 0, then the spreading happens and u(t,x) tends to a positive T-periodic 
function for any moving parameter ^ > 0 regardless of sizes of the initial data uo{x) and self¬ 
limitation coefficient b{t,x)] 

(ii) if Ai(/io; d, a) > 0, we can find a critical value b) of the moving parameter, such that 

the spreading happens and u{t, x) tends to a positive T-periodic function when ^ > fi*{uo, b), while 
the vanishing occurs and u{t,x) will extinct in the long run for /r < fi*(uo,b). 

From our proof we can see that the number ^*{uo,b) is decreasing in uo{x) and increasing in 
b{t,x). The above conclusion (ii) tells us that, in the case \i{ho]d,a) > 0, whether the species 
establishes itself successfully or not depends on the sizes of the initial data and self-limitation 
coefficient. This seems to match the reality better and is supported by numerous empirical evidences 
introduced above. 

Our conclusions indicate that the above mentioned shortcomings of ()6.ip , (16.31) and (|6.4I) will not 
appear if we use the corresponding free boundary problem to describe the spreading or persistence 
instead of the Cauchy problem or initial-boundary value problem with fixed boundary. Furthermore, 
unlike the Cauchy problem model in which the spreading front is represented by an unspecified level 
set of the solution, the free boundary model gives an exact location of the spreading front x = h{t) 
for any given time. We also remark that the important feature of (16.3p . namely the spreading front 
invades at a linear rate in time, is retained by the free boundary model. 

We believe that the cause of the above mentioned defects for problems (|6.ip . (|6.3p and (16.4p is 
that the speed of the expanding front is infinity (too fast) when we use (|6.1I) and (16.3p . and is zero 
(too slow) when we use ()6.4p . In the model (jl.ip . the speed of the expanding front may be more 
reasonable. 

Acknowledgment: The author would like to thank the anonymous referees for their helpful 
comments and suggestions. 
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